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The mass of the bottom quark (both the pole mass Alt and the MS mass mi,) and the strong coupling constant 
as have been determined M from QCD moment sum rules for the T system. In the pole-mass scheme large 
perturbative corrections resulting from coulombic contributions have been resummed. The results of this analysis 
are: Mb = 4.60 ± 0.02 GeV, mb{mt) = 4.13 ± 0.06 GeV and as{Mz) = 0.119 ±0.008. 

1. INTRODUCTION 



A short-distance description in terms of quarks 
and gluons is well suited for inclusive quantities, 
where no reference to a particular hadronic state 
is needed. The vacuum polarization 

iq^.q. - 9,.q') ^{q') = ijdx e''^^T{j, {x)MO)}} 

induced by the vector current = 67^6 can 
be calculated theoretically within the Operator 
Product Expansion (OPE), whereas its imaginary 
part can be experimentally determined from the 
e^e^ — > bb cross-section: 
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Using a dispersion relation the nth derivative 
of H(s) at s = can be expressed in terms of the 
nth integral moment of R{s): 
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where v = a/1 — AM^/s. Mi, corresponds to the 
pole of the perturbatively renormalized propaga- 
tor, whereas the running quark mass in the MS 
scheme renormalized at a scale /i will be denoted 
by mh(p). 

Under the assumption of quark-hadron duality, 
the moments can either be calculated theo- 
retically in renormalization group improved per- 
turbation theory, including non-perturbative con- 
densate contributions, or can be obtained from 
experiment. In this way, hadronic quantities like 
masses and decay widths get related to the QCD 
parameters as, mf, and condensates. 

For large values of n, the moments become 
dominated by the threshold region. Therefore, 
they are very sensitive to the heavy quark mass. 
Owing to the large size of the (ag-^/n)'' Coulom- 
bic corrections, the large-n moments can also be 
used to get a determination of as from the exist- 
ing data on T resonances |Q . 

2. EXPERIMENTAL INPUT 

The first six Upsilon resonances have been ob- 
served. Their masses are known rather accu- 
rately, and their electronic widths have been mea- 
sured with an accuracy which ranges from 3% for 
the T(IS') to 23% for the T{6S). For our pur- 
poses, the narrow-width approximation provides 
a very good description of these states, because 
the full widths of the first three T resonances 
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are roughly a factor 10"^ smaller than the corre- 
sponding masses, and the higher-resonance con- 
tributions to the moments are suppressed: 
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and |] 52 = 



where Tks = T[r{kS) ~* e+ 
1.07a2. 

The e+e^ — s- 66 cross-section above threshold 
is unfortunately very badly measured [Q. The 
second term in Eq. (g) accounts for the contribu- 
tions to Rb above the sixth resonance and is ap- 
proximated by the perturbative QCD continuum. 
The continuum threshold ^/sq should lie around 
the mass of the next resonance, which has been 
estimated in potential models. We have used 
y/sQ = 11.2 ± 0.2 GeV; the lower value includes 
the mass of the sixth resonance and should be 
a conservative estimate. The contribution from 
open B production above the BB threshold and 
below So is very small 0| and has been included 
in the variation of sq. 

The numerical weight of the heavier resonances 
in (||) decreases strongly for increasing values of 
n. The contribution of the T(5S') [T(65)] is 9.5% 
[4%] at n 0; 1% [0.3%] at n = 10; and a 
tiny 0.08% [0.02%] at n = 20. Therefore, taking 
n ^ 10, the uncertainties associated with the con- 
tributions of higher-mass states are very small. 



3. THEORETICAL CALCULATION 

3.1. Perturbation Theory 

The vacuum polarization n(s) can be expanded 
in powers of the strong coupling constant: 

n(s) = n(°)(s) + an(i)(s) + n(2)(s) + . . . , 

with a = as/n. Analogous expansions for R{v) 
and Mn can be written. For the first two terms, 
analytic expressions are available [||j^. 

n(^)(s) is still not fully known analytically. 
However, the method of Pade approximants has 
been recently exploited to calculate H^^^ numeri- 
cally, using available results at high energies, an- 
alytical results for the first seven moments 
{i — 1, . . . , 7) and the known threshold behaviour 



of R^'^\v) [0,11 . This information is good enough 

to obtain an accurate numerical evaluation of 

(2) 

Mil ' for values of n not too large. 

The contributions from diagrams with inter- 
nal quark loops to the spectral function are fully 
known. This allows to check the accuracy of the 
Pade approximation, which for n = 20 is found 
to be better than 10~^. 

The numerical stability of the results can be an- 
alyzed, either using different Pade approximants 
with the full set of information, or constructing 
Pade approximants with one order less by remov- 
ing one datum. For n < 20, the resulting numeri- 
cal uncertainty is below 0.02%, being completely 
negligible for our application. 

The reliability of the Pade approximation has 
been further corroborated in Ref. [^, where the 
first seven terms of the expansion of n^^-* (s) in 
powers of / s have been computed. 

3.2. Coulomb Resummation 

For large n the higher-order perturbative cor- 
rections to the moments grow with respect to the 
leading order. At n = 8 (n = 20) the first order 
correction is roughly 120% (200%) of the lead- 
ing term whereas the second order contribution 
is 140% (340%). This behaviour of the perturba- 
tion series originates from the fact that the rel- 
evant parameter in the Coulomb system is as/v 
which leads to a as\pn dependence of the mo- 
ments [|l0| , [ll|| . Thus, for higher n the perturbative 
corrections become increasingly more important 
and have to be summed up explicitly. 

The resummed spectral function resulting from 
the imaginary part of the Green function for the 
static Coulomb potential is well known [^jj^ . Sub- 
tracting the zero-order contribution, already in- 
cluded in = |i;(3 — w^), it has the form: 



Xy 



(1 



-XV I V 



(3) 



where xy = ir^Cpav and ay is the effective cou- 
pling in the QCD potential; ay is independent of 
the renormalization scale /ia, but it does depend 
on the three-momentum transfer between the 
heavy quark and antiquark, = 4u^M^/(l — w^). 
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The expansion of ay in terms of the MS couphng, 

is known to 0{a^) 

Rc resums the leading (a/w)" and some of the 
sub- leading corrections The corresponding 

terms have to be subtracted from i?'^^ and i?^^^. 
Thus, we can rewrite the perturbative expansion 
of the spectral function in the form 

R{v) = S{a) + Rc + R'^'^a + a^} , (4) 

with 



Sia) 



1 - ACpa + IGCj-a^ 
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(5) 



^(2) ^ ^(2) ^ ^c^j^il) _ 3^ _ 9 ^2^.^ fi) 
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We have factored out in S{a) the correction to the 
vector current "— 4Ci?a", which originates from 
transversal, hard gluons; we have added a term 
IGCpO^ in order not to generate additional cor- 
rections of C(a^) proportional to i?'^*'^. 

After performing the Coulomb resummation, 
the large-moment behaviour of the remaining 
terms is much weaker. For instance, although 
Mn^ /Mn'' increases as ^/n, now Mn'' /Mn^ ~ 
l/\/n. One can easily check that, if ay(g^) is 
evaluated with n/ = 4 light quark flavours, all 
contributions to of 0{l/v), 0(1 nv ) and 

0(1) are canceled in (||), such that R^"^^ vanishes 
in the limit u — > 0. Thus, to O(a^), all thresh- 
old singularities are properly taken into account 
through the Coulomb factor Rc- 

The contributions of 0(f ) which determine the 
constant terms in M^'^ /Mn"^ are not known an- 
alytically. Those terms correspond to the current 
correction from transversal, hard gluons. Thus, 
the splitting of the O(a^w) correction between 
i?^^-' and the global factor S{a) is ambiguous. 

3.3. Power Corrections 

The leading non-perturbative correction is the 
gluon-condensate contribution to the massive 
vector correlator; it is known at the next-to- 
leading order m&. 



The relative growth of A^^^^a/A^n*' is propor- 
tional to n^. Therefore, the non-perturbative con- 
tribution grows much faster than the perturbative 
moments. In addition, in the pole-mass scheme, 
the next-to-leading order correction to A4n,G^ is 
of the same size or larger as the leading term. 
Because the perturbative expansion for the gluon 
condensate cannot be trusted, we shall restrict 
our analysis to the range n < 20 where its contri- 
bution to the bb moments is below 3%. 

4. RESULTS 

To suppress higher resonances as well as power 
corrections, we have restricted the analysis to the 
range n = 8, . . . , 20. Solving the moment sum 
rules for Mb, we can fit Mb to a constant by vary- 
ing Mb and as{Mb). In Figure |^ the resulting 
values for Mb are displayed as a function of n. 
This illustrates that a constant Mb in the range 
8 < n < 20 really produces an excellent fit. 




Figure 1. Mb as a function of n. 



A compilation of all different contributions to 
the errors on Mb and as{Mb) is summarized in 
Table |l|. The dominant uncertainty, due to the 
unknown higher-order perturbative corrections, 
has been estimated in three different ways: 1) 
the size of the known 0{a^) correction (through 
ay to Rc — only the O(a^) contribution has 
been included in the central values — ; 2) the size 
of the O(a^) term i?^^); and 3) the variation un- 
der a change of the scale at which is evalu- 
ated, in the range Mb/2 < /ia < 2Mb. Adding all 
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Table 1 

Error estimates for A/f, and as(Aff,). 





AMb [MeV] 




statistical 


±9.1 


±0.0097 


0{a^) Coulomb 


±7.3 


±0.0222 


scale /j,a 
continuum 

(aGG) 


±0.9 

I 2.1 
±1.3 

_i_2.6 
=■=2.0 

±5.3 


±0.0082 

_i_0.0287 
=■=0.0196 

10.0037 
=■=0.0027 

±0.0036 




±3.1 


±0.0067 


total 


±13.5 


±0.0269 



Table 2 

Error estimates for mh{mi,) and ag{mh). 





Amt [MeV] 


Aas{mb) 


statistical 


±2 


±0.0044 


scale /Zm 


_i_33 
=■=36 


_i_0.0362 
=■=0.0058 


scale 


_i_49 
=■=31 


_i_0.0251 
=■=0.0280 


continuum 


±1 


_i_0.0035 
=■=0.0026 


(aGG) 


±2 


±0.0023 




±3 


±0.0060 


total 


_i_59 
=■=48 


_i_0.0449 
=■=0.0297 



three would double count the error, because the 
uncertainty in an asymptotic series, such as the 
perturbative expansion, is bounded by the size of 
the last known term. For our final results, we 
have chosen to include the errors of varying the 
Coulomb and the O(a^) terms. 

The entry for the gluon condensate results from 
removing the non-perturbative contribution com- 
pletely. 

Adding all errors in quadrature, we finally get: 



Mb = 4.604 ± 0.014 GeV, 
as{Mb) = 0.2197 ±0.0269, 

which implies 

a.(M^) = 0.1184 ±H1]80 ■ 



(6) 
(7) 

(8) 



We have also investigated the same sum rules 
using the MS quark mass. The convergence of 
the perturbative series turns out to be better in 
the MS scheme. However, we have refrained from 
performing a resummation of the Coulomb cor- 
rections because now the velocity v depends on 
the renormalization scale Hm, used to define the 
running mass mb{l^m)- To restrict the C(a^) cor- 
rections to a reasonable size (< 50%), /im should 
lie in the range /x™ = 3.2 ± 0.5 GeV. 

The separate contributions to the theoretical 
errors have been listed in Table The uncer- 
tainty from higher-order corrections is now due 
to the variation of the scales Hm = 3.2 ± 0.5 GeV 
and 2.6 GeV < Ha < 2mf,. The scale fia should 



not be taken lower than roughly 2.6 GeV because 

otherwise the C(a^) correction ^A„ becomes un- 
acceptably large. Adding all errors in quadrature, 
we get: 



nibimb) = 4.13 ±0.06 GeV, 
a,(mb) = 0.2325± ^[oitr^ 

which implies 

a,{Mz) = 0.1196 ±°-°olll ■ 



5. DISCUSSION 



(9) 
(10) 



(11) 



The resulting values for as{Mz) from the pole 
and MS mass schemes turn out to be in very 
good agreement. This is a further indication that 
the uncertainty from unknown higher-order cor- 
rections is under control. In addition, our results 
TOfc(mf,) and Mb for the 6-quark mass satisfy the 
known perturbative relation between the pole and 
MS masses [17 18|, within the errors. 



Combining both determinations of the strong 
coupling constant as, we find 



asiMz) = 0.119 ± 0.008. 



(12) 



We have not averaged the errors of the two de- 
terminations because they are not independent. 
This result is surprisingly close to the current 
world average [0, as{Mz) = 0.1186 ± 0.0036, 
although the error is certainly larger. 

Our results differ from the ones originally 
quoted in Ref. [|j, showing that the errors were 
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grossly underestimated. The main differences 
stem from our inclusion of 0{a^) corrections, 
which are large, and from our more complete nu- 
merical analysis. (In Ref. the analysis was 
performed expanding in powers of 1 /n and keep- 
ing only the leading term; this is a quite bad nu- 
merical approximation; moreover, the actual ex- 
pansion parameter turns out to be 1/y/n). 

The bottom quark mass values obtained by 
us are in good agreement to previous determi- 
nations from QCD sum rules p0|-p3t and a very 
recent calculation from lattice QCD Evolv- 
ing mhirrib) to the Z peak, it also agrees with the 
DELPHI mb{Mz) determination ||5|. Owing to 
the big sensitivity of the moment sum rules for 
the T system to the quark mass, and the good 
control over higher-order as corrections, our re- 
sult is more precise. This result is lower than the 
one obtained in potential models Q . 

The main uncertainty of our results originates 
in the perturbative series. We have estimated the 
effect of unknown higher-order corrections in a 
quite conservative way. However, it is known that 
the pole mass suffers from a renormalon ambigu- 
ity which has been argued to be of C'(lOOMeV) 
27-29j. Throughout our analysis, the pole mass 
has been defined as the pole of the perturbatively 
renormalized quark propagator. Our determina- 
tion (H) might therefore be subject to additional 
uncertainties which go beyond perturbation the- 
ory but which we cannot assess in a precise way. 

It has been pointed out recently that ad- 
ditional Coulombic singularities, not included in 
Rc, could show up at 0{a?). If true, that could 
generate sizeable corrections to the threshold re- 
gion {v 0) and, therefore, to the very large-n 
moments. However, our results have been ob- 
tained at relatively low values of n, and they are 
very stable in the whole range 8 < n < 20. Thus, 
this kind of higher-order corrections seems to be 
safely included in our quoted uncertainties. 
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